Transient behavior in Single-File Systems 
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We have used Monte-Carlo methods and analytical techniques to investigate the influence of 
the characteristics, such as pipe length, diffusion, adsorption, desorption and reaction rates on the 
transient properties of Single-File Systems. The transient or the relaxation regime is the period in 
which the system is evolving to equilibrium. We have studied the system when all the sites are 
reactive and when only some of them are reactive. Comparisons between Mean-Field predictions, 
Cluster Approximation predictions, and Monte Carlo simulations for the relaxation time of the 
system are shown. We outline the cases where Mean-Field analysis gives good results compared to 
Dynamic Monte-Carlo results. For some specific cases we can analytically derive the relaxation time. 
Occupancy profiles for different distribution of the sites both for Mean-Field and simulations are 
compared. Different results for slow and fast reaction systems and different distribution of reactive 
sites are discussed. 

PACS numbers: 02.70.Uu, 02.60.-x, 05.50.-|-q, 07.05.Tp 



I. INTRODUCTION 

Although systems in nature evolve by obeying phys- 
ical laws, it is in most cases difficult or not feasible 
to describe the system properties accurately since de- 
tails of the microscopic dynamics are not fully known. 
Therefore we usually deal with simplified models for 
these systems of which stochastic models are an exam- 
ple. They are thus described by a reduced set of dy- 
namic variables. Although many exact solutions have 
been found, ||, |3[ ^, ^, |^ the vast majority of stochastic 
models cannot be solved exactly. Many results for equi- 
librium systems 0, ||, ||, |lO[ |ll| have been classified. In 
nature, however, equilibrium is rather an exception than 
a rule. In most cases the temporal evolution starts from 
an initial state which is far away from equilibrium. The 
relaxation of such a system towards its stationary state 
depends on the specific dynamical properties and cannot 
be described within the framework of equilibrium statis- 
tical mechanics. Instead it is necessary to set up a model 
for the microscopic dynamics of the system. Assuming 
certain transition probabilities, the time-dependent prob- 
ability distribution Pt (s) to find the system in configura- 
tion s has to be derived from the Master Equation (ME). 
Solving the ME is usually a difficult task. Therefore, the 
theoretical understanding of non-equilibrium processes is 
still at its beginning. Better understanding of these phe- 
nomena would be an important step as non-equilibrium 
systems exhibit a richer behavior than equilibrium sys- 
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tems. |T||T|, |T4[ ^ |6) 

We investigate non-equilibrium processes for a Single- 
File System (SFS) with conversion. In we have 
already elaborated on the special properties of porous 
structures such as zeolites. The one-dimensional nature 
of the zeolite channel leads to extraordinary effects on the 
kinetic properties of these materials. These structures 
are modelled by one-dimensional systems called Single- 
File Systems where particles are not able to pass each 
other. In jT^ we have focused on the steady-state prop- 
erties of a SFS with conversion. The process of diffusion 
in SFS has different characteristics as ordinary diffusion, 
which affects the nature of both transport and conversion 
by chemical reactions. We are investigating the kinetic 
properties of this system, and, more precisely, we are in- 
terested in the properties of the system before reaching 
equilibrium (the transient or relaxation regime). 

Different methods and techniques have been described 
in the literature to solve the ME exactly. ||, 
In spite of the remarkable progress in the field of 
exactly solvable non-equilibrium processes, the major- 
ity of reaction-diffusion models cannot be solved ex- 
actly. Ijl^, It is therefore necessary to use approxima- 
tion techniques in order to describe their essential prop- 
erties (e.g. Mean-Field Approximation, Cluster Approx- 
imation). pi[| Also, as was already realized by Smolu- 
chowski, | |22| fluctuations and correlations may be ex- 
tremely important in low-dimensional systems where the 
diffusive mixing is not strong. Therefore, these approxi- 
mation techniques can give results that deviate strongly 
from the system behavior. Dynamic Monte Carlo meth- 
ods are used to simulate the system according to the ME. 

Few researchers have concentrated on the properties 
of the system in the transient regime and only studied 
the reactivity of the system in this regime. Moreover, 
few research has been done for an open system where 
adsorption/desorption is present at the marginal sites. 
The reason that many of the analytical approaches fail is 
because of the asymmetry of the system. 

In the present work we focus on the non-equilibrium 



phase properties of SFS with conversion. We study the 
relaxation time of the system (time evolution of the sys- 
tem properties, starting with no particles) for different 
sets of kinetic parameters, lengths of the pipe and dis- 
tributions of the reactive sites. In the transient regime 
we observe that MF results are close to the DMC results 
both for slow and fast reaction systems. We outline the 
cases for which the differences are significant. We com- 
pare with the steady-state situation in which the MF 
was not describing the Single-File effects properly, fl^ 
We analyze the situations when analytical results can be 
derived, and we compare these results with Mean-Field 
(MF) and Dynamic Monte Carlo (DMC) results. We 
look at the relaxation of the total loading, loading with 
different components, occupancies of individual sites for 
various parameters and different distributions of the re- 
active sites. As MF is a coarse approximation, for the 
analysis of profile occupancies we introduce a better ap- 
proximation (Cluster Approximation). We analyse the 
results using different analytic methods such as pair and 
MF approximation. Pair Approximation and MF Ap- 
proximation tend to give similar results due to the high- 
order correlations in the system. We investigate the effect 
of various model assumptions made about diffusion, ad- 
sorption, desorption and reaction on the kinetic behavior 
of the system. 

In section II we specify our mathematical model to- 
gether with the theoretical background for analytical and 
simulation results. In section III A we present the vari- 



ous results for transient regimes for the simplified model 
without conversion. We solve numerically and analyti- 
cally the Master Equation in orde r to g et the relaxation 
time of the system. In section 
ory to simplify the rate equations 



[II B 



we use MF thc- 



J_7| of the system for 
the case when all the sites have the same activity to- 
wards conversion. We present the results obtained using 
DMC simulations for the case with conversion when all 



the sites are reactive in section [II C I , and wh en only 
some of the sites are reactive in section [II C 2. For all 



of these cases we compare the DMC results with MF and 
Pair Approximation results. The influence of the posi- 
tion of the reactive sites is also outlined. The last section 
summarizes our main conclusions. 



all spatial correlations in the system, we introduce an- 
other approximation called the Cluster Approximation. 
We have also simulated the system governed by the ME 
using DMC simulations. 



A. The Model 

We model a Single-File System by a one-dimensional 
array of sites, each possibly occupied by an adsorbate. 
This is the model of diffusion and reaction in a one- 
dimensional arrangements of particles with hard-core in- 
teraction. The sites are numbered 1,2,... ,S. A particle 
can only move to the right or to the left if an adjacent 
site is vacant. The sites could be reactive and unreactive 
and we note with iVprot the number of reactive sites. A 
reactive site is the only place where a conversion may 
take place. 

We consider two types of adsorbates, A and B, in 
our model and we denote with Y the site occupation 
of a site, Y=(*, A, B), which stands for a vacant 
site, a site occupied by A, or a site occupied by a B, 
respectively. We restrict ourselves to the following 
mono- and bi-molecular transitions. 

a) Adsorption and desorption 

Adsorption and desorption take place only at the 
two marginal sites i.e., the left and rightmost sites at 
the ends of the system. 



A (gas) 

Am 

Bm 



+ *m, 

■ A (gas) 
B(gas) 



A, 



where subscripts m denotes a marginal site. Note that 
there is no B adsorption. B's can only be formed by 
conversion. 

b) Diffusion 

In the pipe, particles are allowed to diffuse via hopping 
to vacant nearest neighbor sites. 



II. THEORY 

In this section we give the theoretical background for 
our analytical and simulation results. First we specify 
our model and then we show that the defined system 
obeys a Master Equation (ME). We derive then a set of 
exact rate equations from the Master Equation of the 
system describing the reaction kinetics. We look at the 
properties of the system in the transient regime by solv- 
ing these rate equations. In order to do this we have to 
use approximation techniques. We use MF analysis and 
we derive a set of equations that we can solve numeri- 
cally. Because MF is a strong simplification and neglects 



A„ + *„+J < > *„ + An+l 

B„ + *„+J < > *„ + Bn+l 

where the subscripts are site indices: n— 1, 2, . . . , S-1. 
c) Conversion 

An A can transform into a B at a reactive site. 

In the initial state of the system all the sites are va- 
cant (no particles in the pipe) as we are interested in the 
behavior of the system towards equilibrium. 



desorption 




B. Master Equation 



From this follows the rate equation 



Reaction kinetics is described by a stochastic process. 
Every reaction has a microscopic rate constant associated 
with it that is the probability per unit time that the 
reaction occurs. Stochastic models of physical systems 
can be described by a Master Equation. Q 

By a, /3, we will indicate a particular configuration of 
the system i.e., a particular way to distribute adsorbates 
over all the sites. Pa{t) will indicate the probability of 
finding the system in configuration a at time t and Waf3 is 
the rate constant of the reaction changing configuration 
/3 to configuration a. 

The probability of the system being in configuration 
a at time t + dt can be expressed as the sum of two 
terms. The first term is the probability to find the sys- 
tem already in configuration a at time t multiplied by 
the probability to stay in this configuration during dt. 
The second term is the probability to find the system in 
some other configuration /3 at time t multiplied by the 
probability to go from /? to a during dt. 



P„(t + dt) = (1 - dtY,Wpo.)Pc.{t) + dtY,W^pPp{t) 

(1) 

By taking the limit dt ~^ this equation reduces to a 
Master Equation: 



dPc^jt) 
dt 



Y,[w^Mt)-Wf3Mt)]. (2) 




Analytical results can be derived as follows. The value 
of a property X is a weighted average over the values Xa 
which is the value of X in configuration a: 



(X) =^P„X„ 



(3) 



d{X) 
dt 



a/3 



(4) 



Y,W^pPf,{X^ - Xf,). 



C. Analytical methods 

1. Rate equations 

Starting from the Master Equation (2) and using ex- 
pression (4) the rate equations of the system are de- 
rived. Q We denote by H^ds, W^es, W^is the rate 
constants of adsorption, desorption and diffusion respec- 
tively. For simplicity we assume that the rate constants 
of A and B desorption are equal, and also the rate con- 
stants of A and B diffusion are equal. We denote by (Y„) 
the probability that a particle of type Y is on site n, and 
with (Y„Y„_|_i) the probability that a particle of type 
Y is at site n and one at site n + 1, where Y=(*,A,B). 
The coefficients A„, where n=l, 2,..., S, are 1 if site 
n is reactive and otherwise. The rate equations for a 
non-marginal site are 



^^^^ = Wdiff[-(A„*„+i) - (*„_iA„) + (A„_i*„} 
+ (*„A„+i)] - A„W,.x(A„). 



For (B„) we get similarly 
d(B^ ^ 

(*„B„+i)] A„W„(A„) 



^ VFdiff[~(B„*„+i) - (*„_iB„) -I- (B„_i*„) 



The marginal sites also have adsorption and desorption. 
They can be dealt with similary as the conversion. The 



rate equations for A are 
d(Ai) 



which we approximate. This leads to a so-called Cluster 
Approximation. The closure relation 



dt 

d(As} 
dt 



= M^diff h(Ai*2) + (^lAa)] + W^ads(*i) 
-W^dcs(Ai)-AiW,,(Ai), 

= Wdiff [-(*5-lAs) + (*sA5-l)] + W^ads(*5) 



{XY) = ^ {XYZ) 



(11) 



and the rate equations for B 
d(Bi) 



dt 
d{Bs) 
dt 



should now hold for any approximation. We denote here 
with X the occupation of site n, with Y the occupation 
of site n + 1, and with Z the occupation of site n + 2. 
(7) There are various decoupling scheme used in the lit- 

erature j2[ ^, ||, |l^ as approximations for the n-site 
probabilities. For many of these decoupling schemes the 
closure relation (11) no longer holds. For simplicity we 
*iB2)] — Wdcs(Bi) + AiVFi.x(Ai),will use the simplest Cluster Approximation, called Pair 

Approximation for which only the correlations between 
pairs of nearest neighbors (NN) are considered. 

The decoupling scheme for our pair Approximation is 



= W^diffh(Bi*2) 

= VKdiff[-(*5-lAs) + (As-l*s)] - W^des(Bi) 

+ AsW,^{As). 



(8) 



Note that these coupled sets of differential equations 
are exact, but not closed. 



2. Mean-Field 

To solve this coupled set of differential equations, wc 
need to make an approximation for the two-site prob- 
abilities such as (A„*„+i),(B„*„+i), etc. The closure 
relation 



{X)=J2{XY) 



(9) 



(XYZ) = 



{XY){YZ) 
{Y) ■ 



(12) 



It is straightforward to see that the closure relation (11) 
holds. 



1 



{XY){Y) = {XY). 

(13) 



After decoupling, this system of coupled sets of differen- 
tial equations, consisting of the rate equations for one- 
site and two-site probabilities, becomes closed and can 
be solved numerically. 



should hold for any approximation for these two-site 
probabilities. We denote with X the occupation of site 
n and with Y the occupation of site n + 1. The simplest 
approximation is 



{XY) = {X){Y), 



(10) 



i.e. neighboring sites are considered independent. 

The two-site probabilities then become 
(A„*„+i) = (A„)(*„+i), (B„*„+i) = (B„)(*„+i). @ 
This approximation is called the Mean-Field Approx- 
imation and gives us a coupled set of differential 
equations that we can solve numerically. 



3. Cluster Approximation 

The MF Approximation is a strong simplification be- 
cause it neglects all spatial correlations in the system. 
Because the system we analyze is one-dimensional, the 
correlations might be significant and important. The ob- 
vious possibility to eliminate the weakness of the MF 
approach is to introduce another approximation. 

Instead of using MF Approximation for the two-site 
probabilities, we write down their rate equations (see ap- 
pendix). These equations have three-site probabilities 



D. Dynamic Monte Carlo 

We have seen that we can derive approximate analyt- 
ical solutions to the Master Equation. DMC methods 
allow us to simulate the system governed by the Mas- 
ter Equation over time. We simplify the notation of 
the Master Equation by defining a matrix W contain- 
ing the rate constants W^^, and a diagonal matrix R by 
Raf3 = ^7/3 1 if a = /3, and otherwise. If we put the 
probabilities of the configurations Pa in a vector P, we 
can write the Master Equation as 



dP 

'dt 



(R W)P. 



(14) 



where R and W are time independent. We also introduce 
a new matrix Q, Q(t) = exp[— Rt]. This matrix is time 
dependent by definition and we can rewrite the Master 
Equation in the integral form 



P{t) = Q{t)P{0)+ f dt'Cl{t-t')WP{t'). (15) 
Jo 

By substitution we get from the right-hand-side for P{t') 



A. No conversion 



Pit) = m) 

+ [ dt'q{t-t')WQ{t') 
Jo 

+ [ dt' [ dt"Q{t~t')WQ(t' -t")WQ(t") 
Jo Jo 



]P(0). 



(16) 



Suppose at i = the system is in configuration a 
with probability Pq(0). The probability that, at time 
t, the system is still in configuration a is given by 
Qaa{t)Pa{0) — exp{— Raat)Pa{0) ■ This shows that the 
first term represents the contribution to the probabilities 
when no reaction takes place up to time t. The ma- 
trix W determines how the probabilities change when 
a reaction takes place. The second term represents the 
contribution to the probabilities when no reaction takes 
place between times and t' , some reaction takes place at 
time t' , and then no reaction takes place between t' and t. 
The subsequent terms represent contributions when two, 
three, four, etc. reactions take place. The idea of the 
DMC method is not to compute probabilities Pa{t) ex- 
plicitly, but to start with some particular configuration, 
representative for the initial state of the experiment one 
wants to simulate, and then generate a sequence of other 
configurations with the correct probability. The method 
generates a time t' when the first reaction occurs accord- 
ing to the probability distribution 1 — exp[^Raat]. At 
time t' a reaction takes place such that a new configura- 
tion a' is generated by picking it out of all possible new 
configurations p with a probability proportional to Wa'a- 
At this point we can proceed by repeating the previous 
steps, drawing again a time for a new reaction and a new 
configuration, pa, Eq] 



III. 



RESULTS AND DISCUSSION 



In | |l7[ various results for the system with conversion 
(Wrx 7^ 0) and without conversion(Wrx — 0) were re- 
ported. In case Wrx = we have only A particles in the 
system. The total loading [Q) of the system is defined as 
the average number of particles per site. In case Wrx 7^ 
we have B's as well in the system. In this case, the total 
loading (Q), is the sum of the loading with A particles 
(Qa) and loading with B particles (Qb) 



s s 



n=l 



(17) 



71=1 



Note that the total loading of the pipe for the model with 
conversion is the same as for the model without conver- 
sion. We study the relaxation time of the system 
without conversion and of the system with conversion. 



We are interested in the relaxation time of the system 
(transients). We start with the evolution of the total 
loading over time starting from a system with no particles 
at all. As the total loading is the same for the case with 
and without conversion, we will consider for simplicity 
the case with no conversion first. 

As we can derive a finite set of exact rate equa- 
tions(5,6,7,8) it's not necessary to work with the Master 
Equation in this case. With (X„) the probability that 
site n is occupied we have 

= Wdiff[-(X„*„+i) - (*„_iX„) ^^g^ 

+ (X„_i*„) + (=i=„X„+i)] 

when n is not a marginal site. The two-site probabilities 
can be eliminated by using closure relations, 

(X„X„+i) + (X„*„+i) = (X„), (19) 

that hold in this specific case. The probabilities with 
particles on both neighboring sites cancel and the result 
is 

= Wdiff [(X„-i) - 2(X„) + (X„+i)] . (20) 
For the marginal sites we get 

= Wads [1 - (Xi)] - Wdes(Xi) 
+ Wdiff [(X2)-(Xi)] 

and (21) 



rf(Xs) 



dt 



Wads [1 - (Xs)] - Wdcs(Xs) 

+ Wdiff [(Xs-i) - (X5)] . 



These equations are used for the derivation of the re- 
laxation time. The rate equations are going to be sim- 
plified to a point where they are a set of homogeneous 
linear ODE's with only one parameter apart from the 
system size. Dividing by Wads + Wdcs and introducing 
the dimensionless parameters 



A 



w 



ads 



Wads + Wde 
Wdiff 



(22) 



Wads + Wdos' 
T = (Wads + Wdos)t, 



we get 
rf(Xi) 

dT 

rf(x„) 

dT 

djXs) 
dT 



= (Xeq)-(Xl)+A[(X2)-(Xi)], 

= A[(X„_i)-2(X„) + (X„+i)], (23) 

= (Xeq)-(X5)+A[(X5-l)-(Xs)]. 



We can write these equations in matrix- vector notation 



as 



d(X> 
dr 



-M(X)+v, 



(24) 



where (X) is a vector containing the occupancy proba- 
bihties, M is the matrix of coefficients having the form 



/1 + A -A 
-A 2A -A 
-A 2A -A 



-A 2A 



V 



2A -A 
-A 2A -A 
-A 2A -A 
-A l + A/ 



(25) 



and V is the vector that makes the system non- 
homogeneous. The parameter (Xeq) is the equihbrium 

coverage. 

Finally we can make the set of equations homogeneous 
by working with probabilities for vacancies: i.e., with 



we get 



■1 = -^ 



(26) 



(27) 



In order to solve equation (27), we try the substitution 

y = ae-'^^ (28) 

Taking out the exponential leaves us with 

Ma = Loa. (29) 

Removing the time dependence yields relaxation times. 
We see that we have obtained an eigenvalue equation. 
The eigenvalues uj are the reciprocal of relaxation times 
for the corresponding eigenvectors. The relaxation time 
of the system as a whole(trei) is the reciprocal of the 
smallest eigenvalue. We can get this time by simply nu- 
merically solving the eigenvalue equation for given S and 
A. 

Solving the eigenvalue equation analytically. 

For some special cases analytical expressions for the 
eigenvalues can be given. We consider the ansatz 

an = (30) 



If wc substitute a„ from expression (29) into the equation 
for n not a marginal site, we get 



V2-Z 

z 



This equation has two solutions. 



= (jO. 



(31) 



z± 



4/ 



where 



f = -- 

A 



(32) 



(33) 



There are two cases to be distinguished. If / > 4 
then both solutions are real. Because this is the same 
as cj > 4A, wc will have this for slow diffusion and for 
eigenvectors with fast relaxation (large iv). As we are 
interested in the slowest relaxation (small lj) we will look 
at the other case / < 4 or w < 4A. The two solutions 
are then each others complex conjugate. (Note for the 
following that / > 0.) 

Fast diffusion or slow relaxation. 

In this case we can write the solutions as 



z± = re^'^ 



(34) 



where r and ip are both real. In fact the equation for 
z shows that when z is a solution, then so is 1/z. This 
means that r = 1, or 



z± = e=^'^. 



Substitution in the equation for z yields 
cos (fi = 1 — -f. 



(35) 



(36) 



This has indeed only solutions for < / < 4. 

Because there are two solutions for z. the solution for 
an is a linear combination of these two solutions: i.e., 



Or, 



:cie'"'^ + C2e- 



(37) 



From equation (28) we remark that o„ should be al- 
ways real. This means that ci is the complex conjugate of 
C2 . The coefficients will follow from the equations for the 
marginal sites. There arc two of these cqiiations. The 
equation above for costp is a third equation. We have 
four unknowns (ci, C2, w, and </?), but, as only the ratio 
between the coefficients can be determined, we should 
effectively be able to determine all of them. 

Substitution of the expression for o„ in the equations 
for the marginal sites taking into account that c\ and C2 
are complex conjugate (ci = ci^ + icu), leads to 



2(1 - w + A) cos(</3) - 2A cos(2v?) -2 sin(.^) (1 - w + A) + 2A sin(2v5) 

2cos{Sip){l ~Lu + X)- 2Acos((S' - l)ip) -2(1 - lu + A)sin(S'^) + 2Asin((S' - l)(p) 



ClR 

cii 



(38) 



value of trei for infinitely-fast diffusion from the equation 




-1 -0.5 0.5 1 

ln(A) 



1.5 2 2.5 



FIG. 2: The general solution for ln(l/[ij) as a function of 
ln(A), for (fi G [0, tt), and S = 30. 



This equation only has non-trivial solutions (the trivial 
solution is ciji = cij = 0) when the determinant of the 
matrix equals zero. This leads to the following equation 



8A {-LU + X + 1) sin((S' - 2)ip) 
-4:{-uj + X + lfsm{{S-l)ip) 
+4A2 sin((5 - 3)(^) = 0. 



(39) 



We can eliminate cu using cos tf = 1 — //2. 

Using (fi as parameter in equation (39), cp G [0,7r), we 
can get the A's. Equations (35) and (32) gives us the w's. 
We get in this way a; as a function of A, coupled by the 
parameter (p (see figure ^) . 

Solving the eigenvalue equation for the total loading, in 
the case X —>■ oo 

We solve the eigenvalue equation numerically in or- 
der to get the relaxation time of a system. We want 
to describe how the relaxation time of the total loading 
(trei) depends on system parameters like reaction, ad- 
sorption/desorption, diffusion rate constants and system 
size S. In figure ^ we show the influence of diffusion on 
the relaxation time for a system of size S'=30. Note that 
there are two regimes describing the dependence on dif- 
fusion of the relaxation time trei ■ The first regime is for 
slow diffusion, when trei decreases rapidly with increas- 
ing diffusion, and the second for fast diffusion, when trei 
slowly decreases with diffusion to a limiting value. 

Because diffusion is infinitely-fast all the sites have the 
same probability to be occupied and the system is ho- 
mogeneous. We can then analytically derive the limiting 



rf(X„) _ 2 Wads 



dt 



S 



[1 - (Xn)] 



S 



■(X„). (40) 



The first term is the probability of a particle to be ad- 
sorbed at the two open ends, and the second is the prob- 
ability of a particle to be desorbed at the two open ends. 
The probability of a particle to be adsorbed to one end 
equals the adsorption rate constant (Wads) times the 
probability to have there a vacancy (1-^X„)), while the 
probability of a particle to be desorbed equals the desorp- 
tion rate constant (Wdcs) times the probability to have a 
particle ((X„)). 

From the above equation we get the expression for 
(X„), 



(X„) — (Xoq) 



1 — exp 



The relaxation time is 



trei 



-2{Wads+Wde 



s 



2(Wads + Wdes) 



(41) 



(42) 



This limiting value is the same as the one indicated 
by the convergence of the curves in figure ^. The other 
parameter that influence the relaxation time trei is the 
length of the pipe S. For diffusion very fast, trei increases 
linearly with S (see figure ^) . 




FIG. 3: The logarithm of the relaxation time as a function 
of In(S') for the parameters VKads=0.2, VKdcs=0.8 and different 
Wdiff. 



B. Conversion 

In the case with no conversion, we have derived a set 
of exact equations and we have simplified the rate equa- 



tions to homogeneous linear ODE's. Including conversion 
in our model, the two-site probabilities can not be elim- 
inated and an approximation is needed. We use the MF 
Approximation and we get a coupled set of differential 
equations that we can solve numerically. In figure ^, from 
the MF results for the transients, we observe that there 
are two different behaviors determined by conversion. 

For fast reaction systems (see figure ^) , the relaxation 
time of the loading with A's (treiA) and B's {treis) is 
equal to the relaxation time of the total loading (trei)- 
This means that when Q has reached equilibrium, Qa 
and Qb have also reached equilibrium. 

For slow reaction systems (see figure ^), the total 
loading Q relaxes faster to equilibrium than the load- 
ing with A's (Qa) and B's (Qb). The regime between 
Q reaching equilibrium and Qa and Qb reaching equilib- 
rium we call the reaction limited regime. 

We remark that for slow reaction systems in the 
reaction limited regime, the vacancy probability can 
be replaced with the steady-state expression = 
W^dos/(W^ads + W^dos)- The Set of equations for the case 
when all the sites are reactive then becomes 



[(A„+i) + (A„_i) - 2<A„)] 
[(B„+i) + (B„_i) -2(B„)] 
[(A2)-(Ai)] -IF,,(Ai) 



d(A„) 


_ WdiffWdcs 


dt 


W'ads + I^dcs 




- Wr^{An), 


d(B„) 


VKdiffVFdos 


dt 


W^ads + W^dcs 






d(Ai) 


IVdiffWdcs 


dt 


W^ads + W'dcs 




- W^dcs(Ai) + - 


d(Bi) 


VKdiffVFdos 


dt 


W^ads + W'dcs 




-T4^dcs(Bi), 


d{As} 


WdmWdcs 


dt 


Wads + Wdcs 




-Wdcs(As) + - 


d{Bs) 


WdiffWdcs 


dt 


Wads + Wdcs 




- Wdes(Bs). 



WadsWdes 
Wads + Wdcs ' 

■ [(B2) - (Bi)] + W,x(Ai) 

■ [(A5-1) - (As)] - W,.(A5) 

WadsWdes 



Wads + Wdcs 
[(B5-1) - (B5)] +Wx(As) 



(43) 

We can use these approximate MF equations for the 
problem of relaxation of A and B loadings for the case 
when the total loading has already reached steady-state. 
The equations for A's can be written as 



dt 



:M'(A) +v', 



(44) 



is the vector that makes this systems non-homogeneous 
having the non-zero elements for indices 1 and S. We 
substitute (A) — {A)ss + c in the rate equations (47), 
where c is the vector with the deviations of the site occu- 
pancy probabilities with A's from the steady-state. The 
substitution yields 



dc 



= M'[(A),,+c] 

= M'c-FM'(A), 
= M'c, 



(45) 



because M'(A)g^ + v = by definition. 

The equations in c are homogeneous and the matrix 
of coefficients is the same as the one in the original rate 
equations. We can obtain the eigenvalue equation by 
making the substitution 



c e 



(46) 



where (A) is a vector containing the occupacy probabil- 
ities with A's, M' is the matrix of coefficients and v' 



and taking out the exponential from the equations. The 
relaxation time of the A loading(QA) is then the recip- 
rocal of the smallest eigenvalue. We can get this time by 
simply numerically solving the eigenvalue equation. 

Solving the eigenvalue equation for the case with 
conversion 

For fast reaction systems we have seen that using MF, 
the A and B loadings have the same relaxation, which 
is the relaxation of the total loading Q (see figure |^). 
The relaxation of Qa can be thus derived from exact 
equations (13,14). For slow reaction systems we can de- 
rive analytically the relaxation of Qa and Qb for the 
reaction limited regime. In this case we expect the relax- 
ation of Qa and Qb to be determined only by reaction. 
From the MF results we remark that desorption has also 
a strong influence on the transients (see figure ^) . This 
is happening because, in the case desorption is very high, 
the adsorbed particles at the marginal sites will hardly 
diffuse into the pipe, most of them being desorbed imme- 
diately. As a result, few particles will succeed in getting 
to the middle sites and the residence time of the parti- 
cles near the marginal sites will decrease. The loading 
of the pipe with particles will converge than slower to 
steady-state. In figure || we can see that the analytical 
results for the transients obtained solving the eigenvalue 
equation corresponds to the numerical results obtained 
using the MF Approximation for the case all the sites 
are reactive. In figure ^ we have analytical results from 
the eigenvalue equation and MF results for the system 
in reaction limited regime and for different distributions 
of the reactive sites. We observe that for marginal sites 
reactive the results are similar but differ considerably for 
the middle sites reactive and for the homogeneous distri- 
bution of the reactive sites. 

Analytically, solving the eigenvalue equation for a very 
slow reaction system, we find that the relaxation of Qa 




FIG. 4: a) Time dependence of ln(Q*), ln(QA*) and ln{QB*) for fast reaction systems (Wads = 0.8, Wdcs = 0.2, Wdiff=l; 
Wrx=2). b) Time dependence of the ln(Q*), ln(QA*) and \n{QB*) for slow reaction systems (Wads=0.8, Wdes = 0.2, Wdiff=2; 
Wrx=0.01). c) Comparison between time dependence of \n{Q*), ln(QA*), ln(QB*) and the time dependence of the function 
exp( — Wrx *t). For slow reaction systems, the slope does not depend only on Wrx, but also on desorption. We have marked 
the time dependence of the \n{QB*) when Wdcs=0. In a), b) and c) we have marked with * the absolute value of the difference 
between the current value and the steady-state value of the parameter. 




FIG. 5: a) MF results (dashed line) and analytical results obtained solving the eigenvectors in the reaction limited 
regime (continuous line) for time dependence of Qa for the parameters Wads = 0.8, Wdes = 0.2, Wrx = 0.01, Wdifi = 2.0, S = 30 
and all the sites reactive. The upper dashed line represents the total loading (Q) in the system. The total loading Q gives us the 
information when the reaction limited regime starts, b) MF(dashed) and analytical results obtained solving the eigenvectors in 
the reaction limited regim (continuous line) for the relaxation of Qa. We have marked with * the absolute value of the difference 
between the current value and the steady-state value of Qa. 



{trelA 

) as a function of desorption varies with reaction 
for low desorption rates and converges to a limiting value 
for very high rates of desorption (see figure |^a,b) . The 
dependence of A loading on Wdcs has two regimes, the 
first for low desorption rates when A loading strongly 
decreases with desorption, and the second when the A 
loading is converging to a limiting value and the adsorp- 
tion process takes over the system behavior. 

We remark that the analytical results for Qa obtained 
from the eigenvalue equation don't give the MF peak of 
A's particles accumulated in the transient regime for a 
slow reaction system. 



C. Comparisons with simulation results 

We present now the results obtained for the transients 
using DMC methods for different sets of parameters. We 
compare them with the Mean-Field and Pair Approxima- 



tion results. As for very large pipes the computational 
effort is considerable, we study a system of size S* = 30. 
We have considered separately the sets of parameters in 
Table |. 

The sets of parameters from a) to e) are for the cases 
of low loading and from f) to j) for high loading. The 
parameters in the table describe the following situations: 
a) and f ) for very slow reaction and slow diffusion; b) and 
g) for slow reaction and slow diffusion; c) and h) for slow 
reaction and fast diffusion; d) and i) for fast reaction 
and slow diffusion; e) and j) for fast reaction and fast 
diffusion. 

As the MF results for the total loading are exact, there 
are, as expected, no differences between these results and 
those of the simulation for the total loading. We compare 
the MF results with the simulation results for the case 
with conversion in the case all sites are reactive and when 
only some of the sites reactive. 
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FIG. 6: a) Time dependencies for ln(QA*) for a slow reaction system (W^ads=0.8, Wdcs=0.2, Wrx=0.1, Wdiff = 2, S=30) in 
the reaction limited regime when five left and right marginal sites are reactive. b)Time dependencies for \n{QA*) for a slow 
reaction system (W^ads=0.8, W^dcs=0.2, VKrx=0.1, W^diff = 2, S'^SO) in the reaction limited regime when ten middle sites are 
reactive. c)Time dependencies for ln(QA*) for a slow reaction system (Wads=0.8, W^dGs=0.2, Wrx=0.1, W^diff ~ 2, 5=30) in the 
reaction limited regime when a number of ten reactive sites are homogeneously distributed in the system. In a), b) and c) the 
continous line is for the numerical results obtained solving the eigenvalue equation of the system in reaction limited regime, 
and the dashed line is for the MF results. We have marked with * the difference between the current value and the steady-state 
value of the parameter. 



3 



2.5 




0.5 



I ' ' ' ' ' ' ' 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 

FIG. 7: a) The analytical results for the relaxation of Qa as a 
function on Wdcs for different reaction rates. The continuous 
line is for Wix = 0.1 and the dashed line for Wrx = 0.01, when 
VFdes e [0,2]. 



1. All sites reactive 

We first look at the time dependence of the loading 
with A (Qa) and loading with B (Qb)- From the simu- 



case 


Ws.dB Wdes 


VKdiff 


Wrx 


a) 


0.2 


0.8 


0.05 


0.01 


b) 


0.2 


0.8 


0.05 


0.1 


c) 


0.2 


0.8 


2 


0.1 


d) 


0.2 


0.8 


1 


2 


e) 


0.2 


0.8 


10 


2 


f) 


0.8 


0.2 


0.05 


0.01 


g) 


0.8 


0.2 


0.05 


0.1 


h) 


0.8 


0.2 


2 


0.1 


i) 


0.8 


0.2 


1 


2 


j) 


0.8 


0.2 


10 


2 



TABLE I: Table containing the sets of parameters used for 
the simulations 



lation results in figures ^ and |9| we see several regimes 
for the transients. In the case diffusion is slow, the re- 
laxation time is determined by diffusion. When diffusion 
is fast and the reaction slow, then the relaxation time is 
determined by reaction and when both are fast, relax- 
ation time is determined by adsorption/desorption. For 
all these cases the simulations results for the transients 
match the MF results, except when we have low reaction 
rates and fast diffusion for both low and high loading. 

When diffusion is fast and reaction is slow, MF overes- 
timates the amount of A's in the pipe both for transients 
and for steady-state. DMC and MF results indicate an 
overshoot for Qa both for high and low loadings in the 
transient regime. The overshooting appears as a conse- 
quence of the difference between diffusion and reaction 
rates constants. Because the reaction is slow and diffu- 
sion fast, many A's start accumulating into the system 
and they are only later converted into B's. The moment 
tmax when the peak appears is determined by the ratio 
between Wdis and S, but always shortly after the initial 
moment to and it lasts only a short time. The height of 
the peak depends on the total loading (W^ads/Wdos) (see 
figure and on the ratio between Wrx and Wdiff(see 
figure ^ and ^2|). For the case of slow reaction and fast 
diffusion(f,g,h) in table II, the higher the peak, the lower 
the ratio between Wrx and VFdiff- In the table II we give 
the relative height of the peak (AH/Qa) for different 
and Wdiff , at high loading (Wads = 0.8, Wdcs = 0.2), 
where AH is the height of the peak. 

MF overestimates the height of this peak comparing 
with simulation results, but both DMC and MF results 
converge in the same way to steady-state. In figure |l^ 
we show that MF results corresponds to DMC results for 
transients for slow reaction systems. When difussion is 
slow and reaction is fast, MF predicts very well the DMC 
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FIG. 8: DMC and MF results for the time dependencies of Qa and Qb for the cases a, b, c, f, h in TABLE I when all the 
sites axe reactive. The straight lines visible for cases c and h (fast diffusion-slow reaction) are the MF lines. In the other cases 
the MF and the DMC results are indistinguishable. 
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FIG. 9: Time dependencies for ln(QA*) and ln((5B*) for the cases a and f (slow diffusion-very slow reaction at low and high 
loadings), b and g (slow diffusion-slow reaction at low and high loadings), c and h(fast diffusion-slow reaction at high and low 
loadings), d and i (fast diffusion- very fast reaction at low and high loadings), e and j(very fast diffusion-fast reaction at low 
and high loadings) in TABLE I using MF. In all the figures we have marked with * the absolute value of the difference between 
the current value and the steady-state value of the parameter. 
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FIG. 10: MF results for the time dependence of Qa for 
Q=0.2, 0.5, 0.8, H/rx=0.01. 
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FIG. 11: MF results for the time dependence of Qa for 
Vyrx^O.Ol, 0.1, 1, when Wads=QA and Wdos=0.2. 



results. We have also in this case an overshoot for Qa in 
the transient regime. When we have slow diffusion and 
very slow reaction, the height of the peak increases with 
both Wrx and Wdis (see table II for d, e, f). Also in this 
case the MF corresponds to the DMC results. 

As MF ignores the spatial correlations between NN 
sites and MF gives qualitatively good results comparing 
with DMC results, we conclude that MF is a good enough 
approximation for the case when all the sites are reactive. 
This is confirmed also by the comparison between MF 
and Pair Approximation. The Pair Approximation gives 
the same results as MF. 



2. Only some of the sites reactive 

We want to see how the distribution of the reactive 
sites influence the relaxation of A and B loadings. We 
compare the MF and the DMC results for fast reaction 
and for slow reaction systems, for different distributions 
of the reactive sites. We distinguish homogeneous distri- 
bution of the reactive sites, marginal sites reactive, and 
middle sites reactive). 



FIG. 12 

vydiff=o.o5, 0.1 

0.01. 



MF results for the time dependence of Qa for 
1, 2, 10, when W^i,=0.8, VFdos = 0.2, Wr^ = 



case 




WMff 


AH/Qa 








MF 


Sim 


a) 


0.01 


2.0 


1.5905 


2.9521 


b) 


0.1 


2.0 


2.1590 


2.693 


c) 


1.0 


2.0 


2.9911 


2.9984 


d) 


0.01 


0.05 


1.4365 


1.2761 


e) 


0.01 


0.1 


1.5070 


1.4546 


f) 


0.01 


1.0 


1.6458 


2.654 


g) 


0.01 


2.0 


1.5904 


2.9384 


h) 


0.01 


10.0 


1.4095 


3.141 


i) 


2.0 


2.0 


3.329 


3.305 



TABLE II: Table containing the height of the peak for differ- 
ent Wrx and Wdiff 



We have previously seen that for fast reaction systems 
and all sites reactive, the relaxation of Qa and Qb is the 
same as the relaxation of the total loading Q. From 
we know that the results for the total loading, both for 
transients and for steady-state, can be derived analyti- 
cally from exact equations and these results corresponds 
to the DMC results. When only the marginal sites are 
reactive and reaction very fast, DMC and MF results are 
similar with DMC and MF results for the cases when all 
sites are reactive and when it is a homogeneous distribu- 
tion of the reactive sites in the pipe. When the reactive 
sites are situated in the middle of the pipe, the loadings 
Q, QA7and Qb relax slower to equilibrium than in the 
case all the sites are reactive, because it takes more time 
for the A particles to reach the reactive sites. 

For slow reaction systems and all sites reactive, the 
relaxation of Qa and Qb is slower than the relaxation 
of the total loading Q. For this case we can not derive 
Qa from exact equations. The A loading is converging to 
equilibrium simultaneously with B loading. We compare 
thus DMC and MF results for different distributions of 
the reactive sites, for different reaction rates. We will 
analyze two cases, the first for slow reaction {Wrx = 0.1) 
and the second for very slow reaction (Wrx = 0.01). 

In the first case (Wrx = 0.1), the B loading for ho- 
mogenous distribution and middle site reactive reaches 
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FIG. 13: Analytical and simulation results for slow reaction 
systems, WdifF=2, Wrx^O.l, W^ads=0.8, l¥dos=0.2, when all 
the sites are reactive. The simulation (continous line) and 
analytical (dotted line) results for the \n{QA*) as a function 
on time. We have marked with * the difference between the 
current value and the steady-state value of Qa. 



equilibrium faster than in the case of marginal sites re- 
active. This is happening because for marginal sites re- 
active, the B's are formed near the open ends and, con- 
sequently, they can easily desorb and the equilibrium is 
reached later. When the marginal sites are reactive, the 
residence time of the B's is small and the probability to 
find a B on a marginal site is high. The loading with B's 
is increasing more slowly to the steady-state value be- 
cause of the A's that are in the middle of the pipe. The 
same behavior has the A loading converging to equilib- 
rium simultaneously with B loading. When the reactive 
sites are distributed in the middle of the pipe, because of 
blocking, the B's can not reach easily the open ends, the 
probability to find a B on a marginal site is small, and the 
residence time of the B's in the system is large. The load- 
ing with B's is increasing fast to the steady-state value. 
From figure |l^, we see that Qa and Qb for homogeneous 
and middle sites reactive have the same relaxation. 

In the case of very slow reaction {W^^ = 0.01), when 
the marginal sites are reactive, the B loading increases 
faster than in the case of homogeneous and middle sites 
reactive, till a certain time because A particles can reach 
faster the reactive sites. While more B's are formed, the 
B loading is converging slower because the B's can desorb 
relatively fast from the marginal sites, the residence time 
of the formed B's is smaller than in the case when the 
reactive sites are in the middle of the pipe or homoge- 
neously distributed (see figure |lj). Qb converges faster 
in case the middle sites are reactive than in the case of 
homogeneous distribution of the reactive sites. 

We compare the time dependence of the A loading and 
of the B loading using DMC and MF for slow reaction 
systems and different rates of conversion when diffusion 
is fast (see figure We see that differences appear in 
the transient region as well as in the steady-state for all 
the distributions of the reactive sites but very prominent 
for marginal and middle sites reactive. In this case, for 
homogeneous distribution of the reactive sites the differ- 



ences between MF and DMC are small (see figure 16). 

For the other cases (fast reaction-slow diffusion, fast 
reaction-fast diffusion, slow reaction-slow diffusion), MF 
gives good results compared to DMC for all the distri- 
butions of the reactive sites. In figure ^ we compare 
the MF results for the transients for the case when five 
marginal sites are reactive. 

We finally look at the site occupancy of the pipe. The 
MF profiles for marginal sites reactive show in the tran- 
sient regime accumulation of the A's in the middle of the 
pipe for a slow reaction system and fast diffusion. Be- 
cause reaction is slow, the A's can pass without reacting 
to the nonreactive sites. As the loading is increasing, the 
residence time of the particles increases, and the proba- 
bility to find an A in the middle of the pipe is decreasing 
because of blocking (see figure p^ . 

Dynamic Monte Carlo simulations show the same be- 
havior. In figure ^ we see as well accumulation of A par- 
ticles in the middle of the pipe in the transient regime. 
As a consequence, until the equilibrium is reached the 
middle sites have different contribution to the occupancy 
in the pipe depending on their position. In the steady- 
state all the middle sites have the same contribution to 
the occupancy profiles |17|. 



IV. 



SUMMARY 



We have used DMC and analytical techniques to study 
the properties of Single-File Systems in the transient 
regime. 

We have derived exact equations to solve the relaxation 
time of the whole system (trei)- We found that there are 
two regimes describing the dependence on diffusion of the 
relaxation time trei ■ The first regime is for slow diffusion, 
when trei decreases fast with increasing diffusion, and the 
second for fast diffusion, when trei slowly decreases with 
diffusion to a limiting value. We have analytically derived 
this limiting value of trei for infinitely-fast diffusion. 

We have also studied the transients in the case with 
conversion. MF results show that there are two differ- 
ent behaviors determined by conversion. For fast reac- 
tion systems, the relaxation time of the loading with A's 
(treiA) and B's {treis) IS equal to the relaxation time of 
the total loading (trei)- When Q has reached equilib- 
rium, Qa and Qb have also reached equilibrium. For 
slow reaction systems, the total loading Q relaxes faster 
to equilibrium than the loading with A's (Qa) and B's 
(Qb)- The regime between Q reaching equilibrium and 
Qa and Qb reaching equilibrium we call the reaction lim- 
ited regime. In the reaction limited regime, MF shows 
that not only reaction, but also desorption has a strong 
influence influence on the transients. We find that the re- 
laxation of Qa (treiA) as a function of desorption varies 
with reaction for low desorption rates and converges to a 
limiting value for very high rates of desorption. 

DMC results shows several regimes for the transients, 
in the case all the sites are reactive. In the case diffu- 




FIG. 14: a)Time dependencies for ln(QB*) when W^ds = 0.8, Wdes = 0.2, VFdiff = 2,Wtx = 0.1 using MF for marginal 
sites reactive, homogeneous distribution, middle sites reactive, and all sites reactive. b)Time dependencies for \ii.{Q^) when 
W'ads = 0.8, Wdes = 0.2, Wdiff = 2, Wrx = 0.01 using MF for the same distributions of the reactive sites as in a. In a) and b) we 
have marked with * the absolute value of the difference between the current value and the steady-state value of the parameter. 
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FIG. 15: a) DMC simulation of time dependencies for Qa and Qb for VFads=0.8, Wdes=0.2, Wdifi=0.05, lVrx=0.01, when 
blocks of 5 marginal sites are reactive (S'=30). The straight lines corresponds to the MF results, b) DMC simulation of time 
dependencies for Qa and Qb for Wads=0.8, Wdes=0.2, VFdiff=2, M^rx=0.1, when blocks of 5 marginal sites are reactive (5=30). 
The straight lines corresponds to the MF results. 



sion is slow, the relaxation time is determined by dif- 
fusion. When diffusion is fast and the reaction slow, 
then the relaxation time is determined by reaction and 
when both arc fast, relaxation time is determined by ad- 
sorption/desorption. For all these cases the simulations 
results for the transients match the MF results, except 
when wc have low reaction rates and fast diffusion for 
both low and high loading. In this case MF overesti- 
mates the amount of A's in the pipe both for transients 
and for steady-state. DMC and MF results indicate also 
an overshoot for Qa both for high and low loadings in 
the transient regime, that appears as a consequence of 
the difference between diffusion and reaction rates con- 
stants. 

When only some of sites are reactive, for fast reaction 
systems, MF gives good results compared to DMC for 
all the distributions of the reactive sites. When only the 
marginal sites are reactive and reaction very fast, DMC 
and MF results arc similar with DMC and MF results 
for the cases when all sites are reactive and when it is 
a homogeneous distribution of the reactive sites in the 
pipe. When the reactive sites arc situated in the middle 
of the pipe, the loadings Q, Qaj and Qb relax slower to 



equilibrium than in the case all the sites are reactive, 
because it takes more time for the A particles to reach 
the reactive sites. 

For slow reactive systems, differences between DMC 
and MF results appear for transients for different distri- 
butions of the reactive sites when diffusion is fast, but 
are very prominent for marginal and middle sites reac- 
tive. For homogeneous distribution of the reactive sites 
the differences between MF and DMC are small. For 
slow reaction, wc find that the B loading for homogenous 
distribution and middle site reactive reaches equilibrium 
faster than in the case of marginal sites reactive. For very 
slow reaction, Qb increases faster at the beginning than 
in the case of homogeneous and middle sites reactive and 
becomes slower as the B's arc formed. Qb converges also 
faster in the case the middle sites are reactive than in the 
case of homogeneous distribution of the reactive sites. 
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FIG. 16: a) Time dependencies for \ti.{Qa*) for a slow reaction system (VKads=0.8, Wdes=0.2, VKrx=0.1, Wdis = 2, S=30) 
when five left and right marginal sites are reactive, b) Time dependencies for \n{QA*) for a slow reaction system (W^ads=0.8, 
V[^dGs=0.2, Wrx=0.1, WdiB = 2, S=30) when ten middle sites are reactive, c) Time dependencies for \n{QA*) for a slow reaction 
system (Wads=0.8, H^dos=0.2, Wrx=0.1, Wdm ~ 2, 5=30) when ten reactive sites are homogeneously distributed in the system. 
In a), b) and c) the continous line is for the DMC results and the dashed line is for the MF results. We have marked with * 
the difference between the current value and the steady-state value of the parameter. 
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FIG. 17: Time dependencies for ln((5A*) and ln(QB*) for the cases a and f(slow diffusion-very slow reaction at low and high 
loadings), b and g(slow diffusion-slow reaction at low and high loadings), c and h(fast diffusion-slow reaction at low and high 
loadings), d and i(fast diffusion-fast reaction at low and high loadings), e and j(vcry fast diffusion- fast reaction at low and high 
loadings) in TABLE I using MF when 5 of the marginal sites are reactive. In all the figures wo have marked with * the absolute 
value of the difference between the current value and the steady-state value of the parameter. 




FIG. 18: MF profile occupancies ((An) and (Bn)) for the case of 5 marginal sites reactive before steady-state is reached 
for a slow reaction system of length S = 30 and parameters VFads=0.8 and Wdes=0.2, VFdiff=2, Wrx=0.1. The lines at low 
occupancies correspunds to (An) profile occupancies after 200, 150 and 100 time units in this order from the bottom to the 
top. The lines at high occupancies correspunds to (Bn) profile occupancies after 200, 150 and 100 time units in this order from 
the top to the bottom. 




FIG. 19: DMC results for the occupancy profiles((An) and (Bn)) for slow reaction systems (VFads=0.8, Wdes=0.2, VFdiff=2, 
M^rx=0.1) in the case 5 marginal sites are reactive, before steady-state is reached. The lines at low occupancies corresponds to 
(An) profile occupancies after 280, 100 and 20 time units in this order from the bottom to the top. The lines at high occupancies 
correspunds to (Bn) profile occupancies after 280, 100 and 20 time units in this order from the top to the bottom. 



VI. APPENDIX 



A. Derivation of the equations for two-sites 
probabilities 



In order to solve the rate equations of the system we 
use Cluster Approxiniation. For simplicity wc consider 
an approximation that considers only the correlations be- 
tween pairs of NN sites - Pair Approximation. To the al- 
ready derived rate equations for onc-sitc probabilities we 
add the two-site probability equations. We write write 
these equations in terms of three sites probabilities and 
we use then the decoupling scheme 



{XYZ) 



{XY){YZ) 
{Y) ■ 



We give here the equations for the two-site probabili- 
ties for A and * occupancy of these sites. The equations 
for the left marginal sites are 



"^^^^d/^^^ = W^diff *s-i As) - {*s-2As-iAs)) 

- 2Wr^{As-iAs) + W^,{As-i*s) 

- Wd,,{As-iAs) 

"'^^^^"Z*^^ = WM{As-i *s *s+i) + {*s-iAs) 

- {As-i*s) - (*S-2^S-l*s)) 

- Wr^{As-l*s) - Wads(^S-ll*S> 

+ Wd,,{{As-iAs) + {As-iBs)) 



d{*s-iAs) 
dt 



■■ Wdis{{As-i*s) + {*s-2As-iAs) 

+ {*s-2Bs-iAs) - {*s-iAs) 

- {As-2 *s-i As) - {Bs-2 *s-i As)) 

- Wrx(*S-l^s) + W^ads(*S-l*S> 

es {*s-iAs) 

Wdiff((*S-2-4s_i*s) + (*S-2-Bs-l*S> 

- {As-2 *s-i *s) - (-Bs-2 *s-i *s)) 

- Wads(*S-l*s) 

+ Wdes((*5-1^5) + {*S-lBs)), 

(49) 

And, finally, the equations for the non-marginal sites 



(47) d{*s-i*s) 



dt 



d{A^A2) 
dt 

d{Ai*2) 
dt 



VKdiff ((Ai *2 ^3) - (^1^2*3)) - 2W,^{AiA2) 

+ M^ads(*1^2) - Wacs{AiA2) 

- W^diff((*1^2) + (^1^2*3) + (^1^2*3) 

(^1*2) - {Ai *2 ^3) - (^1 *2 S3)) 
W;x(^l*2) + Wads(*l*2) 
Wdo,{Ai*2) 



d{AnAn+l) 

dt 



d{An*n+l) 

dt 



rf(*1^2) 

dt 



VPdiff ((^1*2) + (*1 *2 A3) - {*lA2) - (*lA2*3)) d{*nAn+l) 



rf(*i*2) 
Jt 



- Wrx(*lA2) - Wads(*1^2) 
+ Wa,s{{AiA2) + {B^A2)) 

H^diff((*1^2*3) + (*1-B2*3) 

- (*1 *2 ^3) - (*1 *2 -B3)) 

- W'ads(*l*2) + W^des((A*2) + (-81*2)), 



dt 



(48) 



where {A^ *2 A3) = i^iislgsM , etc. 

Almost similar are the equations for the right marginal 
sites 



Wdiff((A„ *„+i An+2) + {An-l *n An+l) 

- {*n-lAnAn+]) - (A„A„+i*„+2)) 

- 2Wpx(A„^„+i) 

= W^diff((* ) + {An-l *n *n+l) 

+ {AnAn+-L*n+2) + {AnBn+\*n+2) 

- {An*n+l) - {An *n+l An+2) 

{An *n+l Bn+2) - {*n-lAn*n+l)) 

- W^^{An*n+l) 

Wdia{{An*n+l) + (*n An+2) 

+ {*n-lAnAn+l) + {*n-lBnAn+l) 

- {*nAn+l) — {*nAn+l*n+2) 

- {An-l *n An+l) - {Bn-1 *n An+l)) 

- Wr^{-*nAn+l) 

= Wdifi'((*„Ai+l*n+2) + {*nBn+l*n+2) 

+ {*n-lAn*n+l) + (*n-l-B„*„_|_i) 

- (*n An+2) - {*n *n+l Bn+2)) 

- {An-l *n *n+2) - {Bn-1 *n *n+l)) 

(50) 



We have the possibility to determine the one-site proba- 
bilities by calculating the sum of the two-sites probabili- 



d{*n*n 



+11 



dt 



ties 

^{XY) = {X). 

Y 
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